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Abstract 

We investigate the equilibrium of a fluid in contact with a solid boundary through a density- 
functional theory. Depending on the conditions, the fluid can be in one phase, gas or liquid, or 
two phases, while the wall induces an external field acting on the fluid particles. We first examine 
the case of a liquid film in contact with the wall. We construct bifurcation diagrams for the film 
thickness as a function of the chemical potential. At a specific value of the chemical potential, 
two equally stable films, a thin one and a thick one, can coexist. As saturation is approached, the 
thickness of the thick film tends to infinity. This allows the construction of a liquid-gas interface 
that forms a well defined contact angle with the wall. 

Keywords: contact angle; density-functional theory. 
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I. INTRODUCTION 



Wetting phenomena have received considerable attention over the past few decades (com- 

f] fl 

prehensive and detailed reviews are given by |6| and [9|). Central to any description of 
wetting is the presence of a contact line which involves a three-phase conjunction (gas- 
liquid-solid). In the case of dynamic wetting the associated moving contact line problem is 
characterized by a stress singularity at the contact line. As a consequence, formulating this 
problem in the framework of conventional fluid mechanics leads to fundamental difficulties 
concerning the modeling of the contact line region 

ana 

. Accordingly, a variety of models 
have been proposed to alleviate these difficulties and to address the behavior of the associ- 
ated dynamic contact angle, with the most popular approaches being the replacement of the 
no-slip condition with a slip model or the elimination of the contact line all together with 
the use of a thin precursor layer. More recently, new approaches/theories have appeared, 
based e.g. on local diffuse-interface/Cahn-Hilliard-type approaches which, in general, cannot 
be rigourously justified. Other recent approaches include hybrid molecular-hydro dynamic 
simulations. However, the question of how to rigorously coarse grain from the micro- to 
the macro-scale and to accurately transfer molecular information to the macro-scale has not 
been adequately addressed as of yet. 

The difficulties with the contact line region stem from the multi-scale nature of the 
problem. At the macroscale the usual laws of hydrodynamics apply. At the nanoscale 
intermolecular interactions, often described by molecular dynamics/Monte Carlo simula- 
tions, dominate. Yet, phenomena occurring at the nanoscale often manifest themselves at 
the macroscale (where individual molecules are not considered or equivalently the hydrody- 
namic laws consider a very large number of molecules at the same time). Despite drastic 
improvements in computational power, molecular simulations are still only applicable for 
small fluid volumes. 

A compromise between conventional hydrodynamics and molecular simulations can be 
achieved by density-functional theory (DFT). On the one hand, DFT is able to retain the 
microscopic details of a macroscopic system at a computational cost significantly lower than 
that used in molecular simulations. On the other hand, DFT is rigorous compared to con- 
ventional phenomenological approaches. It is applicable for both uniform and non-uniform 
(density exhibits spatial variation) as well as confined systems (e.g. in the presence of a wall) 
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within a self-consistent theoretical framework provided by the statistical mechanics of fluids. 
At the same time, substantial progress has been made in recent years in the development 
of realistic free-energy functionals which take into account the thermodynamic non-ideality 
attributed to the various intermolecular forces. The papers by js], 23 1 and [l| outline a 
DFT framework currently widely accepted by the statistical mechanics of fluids commu- 
nity and which has also been used successfully in describing equilibrium configurations in 
many different contexts, from fluids in pores to liquid crystals, polymers and molecular 
self-assembly. 

Several studies have examined the equilibrium of liquids in contact with solids in different 

settings and in the framework of DFT/statistical mechanics of fluids. For example, 15 1 

□ 

adopted a DFT approach based on a free-energy functional proposed by [12]. This study 
focused primarily on a one-dimensional (ID) configuration, namely a liquid film in contact 
with a planar solid substrate but it did discuss contact lines; for a example it gave an 
approximate expression for the equilibrium contact angle obtained from the sharp-interface 
imit of the DFT approach adopted and for large distances from the solid substrate. Recently, 
3|-|5[ adopted a DFT approach based on the framework described by fl, Q and to 
calculate nanodrops on chemically/physically inhomogeneous inclined/planar substrates and 
to obtain the dependence of the contact angle of nanodrops on planar horizontal substrates 
on the parameters of the intermolecular interactions. 

Our aim here is to provide a rigorous methodology for the treatment of contact lines. 
Our approach is also based on the DFT framework described by Q], 23 1 and [ij]. We focus 
on the equilibrium of a fluid in contact with a planar solid substrate whose understanding 
is essential for the substantially more involved dynamics. We first examine in detail the 
ID case of a liquid film in contact with the substrate. Particular emphasis is given on 
the bifurcation diagrams for the film thickness as a function of the chemical potential. 
This is a necessary step prior to understanding the more involved two-dimensional (2D) 
case, i.e. the three-phase contact line. We subsequently focus on the prewetting transition 
occurring at a specific value of the chemical potential where two equally stable films, a 
thin one and a thick one, coexist. The two thicknesses are connected through a ridge- 
like interface which is sufficiently smooth in the vicinity of the thin film for a long-range 
wall potential but it becomes steep there when the wall potential is a short-range one. 
When the co-existence value of the chemical potential equals to the saturation one, the 
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thickness of the thick film tends to infinity. This is the case of a hquid wedge in contact 
with the sohd substrate and with a well-defined three-phase contact line a few molecular 
diameters from the substrate. The wedge seems to persist for all distances from the substrate. 
Hence, even though DFT is a microscopic approach, it allows for the construction of a 
macroscopic quantity such as contact angle. Moreover, unlike macroscopic approaches e.g. 
Young's equation which naturally require information on macroscopic parameters, such as 
the surface tensions between the different phases, DFT relies only on first principles, namely 
information related to intermolecular parameters. It also elucidates the structure of the 
three-phase contact line region and the precise role of fluid- fluid and fluid-solid (long-/short- 
range) interactions there. 

II. PROBLEM DEFINITION 

A. Setting 

We consider part of a simple fluid inside a volume V in contact with a horizontal planar 
substrate as shown in figure [H The system, whose boundary is denoted with the closed 
dotted line in the figure, is open, i.e. fluid particles can come in and out, and its surroundings 
are at temperature T and chemical potential fi. Depending on the conditions, the fluid can 
be in one phase (gas or liquid) or two phases. The fluid particles interact through a pair 
potential 0(ri2), where is the distance between the centers of mass of the particles. The 
main effect of the substrate is to induce an external field V^(r) acting on the fluid particles 
with r the position vector of the inertial center of fluid particles. For the sake of simplicity, 
gravity is neglected. A Cartesian coordinate system (x, y, z) is chosen such that x and y are 
parallel to the wall surface while z is the outward-pointing coordinate normal to the wall. 
Of particular interest is the region around the solid/fluid interface. 

B. Density-functional theory of fluids 

The equilibrium state of the system is characterized by the fluid density no(r), in units 
of number of particles per unit volume (number density), defined as jsj: 
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FIG. 1: Sketch of the physical setting. 



with r G V. To obtain nQ{r), we use elements from equilibrium statistical mechanics of 
fluids, in particular DFT. It has been shown that there exists a functional Q[n\, defined on a 
set of functions n{r) compatible with the external potential V^(r), which has the property to 
be at a global minimum (for a given temperature, chemical potential and external potential) 
when n(r) is equal to the equilibrium density profile no{r) [s, [l^ (minimum principle). A 
second functional F[n] can also be introduced such that Q[n] reads: 



where the integral is understood to be a volume one over V. The intrinsic free energy and the 
grand potential of the system at equilibrium are then equal to -F[no] and fi[no], respectively 
(such quantities can only be defined at equilibrium). 

By using variational calculus, it can be easily shown that a necessary condition for a 
minimum of Q[n\ is: 



This equation is supplemented by appropriate boundary conditions, i.e. the value of n outside 
V. In our case, the external potential V{r) is due to the interaction between the fluid and 
the wall. Note that since our system is open (grand-canonical ensemble formalism), the 
total number of particles, N, is not given, but it is obtained instead once no(r) is known: 




(2) 




(3) 



N = J druQ^r). For simplicity, we drop the subscript from n hereinafter. 
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C. Fluid modeling 



Several approximations for F[n\ have been proposed over the years. One which has proved 
to be successful in a number of cases is based on a perturbation approach. The basic idea 
is to split the pair potential 0(^12) into two terms: 0(^12) = (prijii) + 0p('^i2)- The first 
term, 0r(^i2), is a reference potential and usually corresponds to the harshly repulsive part 
of 0(ri2), while the second term, </'p(ri2), acts as a "perturbation" to </'r(?"i2) and is in general 
a slowly attractive potential. Such separation, which is reminiscent of the approach behind 
the van der Waals equation of state for gases and can thus be viewed as a generalization 
of this approach, is in general well-suited for treating dense fluids. Its underlying physical 
motivation is that the structure of a fluid is determined mainly by the repulsive forces. After 
some additional simplifications we have the following compact expression for F[n\ jsl. Il6|: 

F[n] = F,[n] + \ jj drdr' n(r)r2(r')0p(|r - r'|). (4a) 

The first term of equation (Hal) accounts for (j)r{ri2) and can be further simplified by using 
the co-called local density approximation: 

F,[n] = Jdrn{r)f{n{r)), (4b) 

where / is the local free energy per particle of the reference fluid. The main simplification 
used to obtain the second term of equation (Hal) , which involves the attractive part of the 
interaction potential 0(ri2), is to express the pairwise distribution function, i.e. the proba- 
bility that two particles occupy the positions ri, r2, as the product of ?7(ri), ?7(r2) and the 
radial distribution function of the uniform reference fluid. 

The decomposition of the pair potential in two parts can be done in different ways. In 
the Barker-Henderson approach for example the attractive part is. 
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ri2j \r12J 



if ri2 < a 

(5) 

if ri2 > (T, 



where e is a parameter that measures the strength of the potential. Another possibility is 



the decomposition suggested by 25|. Here we adopt the Barker- Henderson approach. The 
reference potential is then, cf)^ = cf) — 0p. However, to take full advantage of the perturbation 
technique, the reference system should be a well-known fluid: the preferred choice is often 
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a hard-sphere fluid, whose bulk local free energy is well described by the Carnahan-Starling 
expression: 



/3/(T, n) = HA'n) - 1 + ^^5^^ (6a) 

(1 - r]) 



where 



V = ^d'n, (6b) 





is the packing fraction, d is the hard-sphere diameter, A is the thermal de Broglie length 
and /3 = l/{kBT) with ks the Boltzmann constant. This means that (pr is in fact approxi- 
mated by a hard-sphere potential. The associated molecular diameter d, which appears in 
equation 
potential 

Instead of equations (]4bl) and ( l6al) . one could use a more refined approach for a hard 



6b]) and hence fj6a|l . can be linked to a and the reference part of the interaction 
2| but for simplicity we assume d = a. 



sphere fluid based on Rosenfeld's fundamental measure theory [18|. In general, this approach 
gives better results for the fluid density at small distances from the wall (a few molecular 
diameters) but requires a more involved computational treatment. Our aim here is to keep 
the formalism as simple as possible and yet retain the basic ingredients of the underlying 
physics. 

By substituting equations (Ha|) and (llbll into equation ([2]), we obtain: 

/ir(n) + y"dr'n(r')0p(|r-r'|) + r(r) = /i (7a) 

where 

_ f d{nf{T,n)) \ 
^-^"^=1 dn J,' ^^^^ 
is the chemical potential of the reference system which in the case of a hard sphere fluid and 



/3/Xr(n) = In (A nj + -3 . (8) 

[i-v) 



using equation ( I6al) reads: 



When the external potential vanishes and the fluid is uniform, i.e. its density does not 
exhibit any spatial variation, equation ( Ha|) combined with equation ( l4b|l . reduces to 

F[n] = Vnf{n) (9a) 

where 

-j^ r + oo 

f{n) = f{n) — an and a = — / dr0p(|r|). (9b) 

2 J -00 



Equation ([7]) is then an algebraic equation with parameters T and fi. Using for / the 
expression given in equation f l6a|) . the usual thermodynamics for liquid/gas systems applies. 
The solutions to equation ([7]) correspond to extrema of the functional Q[n\. There is one 
solution for T > and one to three solutions for T < T^, where is the critical temperature 
(/3c = 1/ (AjbTc)). In the latter case, there are two minima of equal depth for Q, at small and 
large densities, respectively, and a maximum at intermediate densities. The middle solution 
is unstable while the other two correspond to the liquid and gas bulk densities, denoted as 
UgasiT,^) and niiq(T, /i), respectively, with n^asiT^jji) < niiq(T, /i). The chemical potential 
at the gas/liquid transition will be denoted as /isat(^)- The gas is the preferred state for 
< fJ'sa.tiT). At fi = fissLtiT) the gas and liquid bulk are at equilibrium while for /i > fisa.t(T) 
the liquid is the preferred state. Note that in all cases considered here T < as close to 
the critical temperature thermal fluctuations, which are not taken into account in the DFT 
formalism adopted here, become significant. 



D. Wall potentials 

Three types of wall potential are considered in this study. All are attractive at large 
distances. The first two have the same infinite repulsive part but they differ in the range of 
the attraction term. The short-range potential reads. 



Vsr{z) 



while for the long-range one we use. 



-oo 



-e„ exp 



-oo 



2/ 2/\f 



if z < 
if z > z^ 



if 2; < 2;„ 
if z > Zw. 



(10a) 



(10b) 



In these expressions, cr^ and z„ are three wall parameters; the first two are strictly 
positive and are related to the strength and range of the potentials. The expressions in 
equations llOal and llObI are useful in that they allow us to examine the effect of the wall 
attraction, and in particular its range, on the fluid equilibrium state. The last wall potential 
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has a smooth repulsive part compared to the previous ones (i.e. it is non-infinite): 

'OO if z < 

2e 



(lOc) 

if z > Zw. 



15 V 

It can be derived by considering that the wall is made of a uniform density of particles 
interacting with the fiuid particles through a Lennard- Jones potential (see Appendix A). By 
setting cr„,LR = ^w,lr - z„,lj and e^.LR = (crw,Lj/crw,LR)^ew,Lj, the two last potentials have 
the same attractive part. Finally, it is also possible to have a purely repulsive wall with no 
attraction (a hard or dry wall), 



Viz) 



+00 if 2; < 2;„ 
ii z > z^. 



III. ID PROBLEM 
A. Equations 

In this section we assume that the system is invariant along the x and y directions and 
infinite in the ^-direction. Equation (!7al) then becomes: 

/+00 
dz'n{z')^id{z' - z) + V{z) = (11a) 
-00 

where 

$id(^) = jj °dxdy0p(v/x2 + y2 + z2j. z = z'-z. (lib) 

The integral in (lllbl) can be performed analytically (see Appendix A). To solve numer- 
ically the integral equation ( ITT]) a simple iterative procedure can be employed based on a 
Picard scheme. The integral is first computed with the values of the density obtained from 
the previous iteration by using a simple trapezoidal rule and then the remaining terms are 
solved for ra, which gives the new value for the density. We found that, with the exception 
of unstable solutions, this scheme is robust and allows the use of an initial guess which can 
be far from the solution. However, convergence tends to be slow after the first iterations. 
This can be greatly improved by using instead a more involved scheme based on a Newton 
procedure in which the Jacobian matrix is appropriately simplified. Details are given in 
Appendix B. 



B. Density profiles 



In figures [2] and [31 we present typical profiles obtained by solving equation (fTTIl . Figure [2] 
depicts the density profile of the fluid along the 2;-direction when /i = fisat{T) in the absence 
of the wall. For this value of fi, both bulk liquid and bulk gas are equally possible states {fl 
has two variational minima of equal "depth" at nuq^gas)- The interface region is just a few 
molecular diameters thick. The location of the sharp interface, say zi, can be obtained from 
the Gibbs dividing surface defined from /f'^(7^iiq — n)dz = f^°°{n — nga,s)dz. 

Figure [3] shows profiles obtained when the wall is switched on and exerts a long-range 
attraction on the fluid molecules. In this case, fi < /isat(^) so that the bulk gas is the more 
stable phase in the absence of the wall {Q has a single variational minimum at rigas). As 
z — 7- +00 the density tends to the density of the gas, rigas, but near the wall the attractive 
potential induces a bump in the density profile. If /i is close enough to fisa.t(T) (solid hne 
in figure [3]) , this bump is more pronounced while the fluid density in that area is roughly 
equal to the metastable liquid bulk density (corresponding to a local variational minimum 
of n). A thin liquid film is effectively formed between the wall and the gas. It can only exist 
due to the attraction with the wall, in other words an attractive wall stabilizes a thin liquid 
film in contact with it. We also note the small oscillations in the density profile near the 
wall corresponding to adsorption of the liquid particles there. More pronounced oscillations 
would have been found in this area if the reference system of hard sphere fluid had been 
treated with a more refined model like the Rosenfeld fundamental measure theory mentioned 

n 

earlier [18|. Also, when a thin film is present, we observe that the interface between the film 
and the gas phase is very similar (although /i is different) to the profile given in figure [2] 
suggesting that in that case, the wall has little influence. In fact, as /i — /is^t; the thickness 
of the liquid film tends to infinity and we approach the case depicted in figure [21 i.e. it is 
like the wall is not even present (with the exception of course of the area close to it where 
the density oscillations occur). 

For given external conditions of temperature and chemical potential, an illustration of 
the dependance of the density profile on the characteristics of the wall potential is shown 
in figure [H We first note when comparing wall potentials Vsr (dashed-dotted line) and 
Vlr (dashed line), that the liquid film thickness in the first case is considerably smaller 
even though all wall parameters are identical and in particular the integrals between z„ and 
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FIG. 2: Density in the vicinity of a liquid-gas interface for = /igat at T = 0.7 Tc in the absence 
of the wah. The domain size used in the computation is [0-50] and the number of points is 1000. 
The bulk liquid density is nnqd^ ~ 0.684 and the bulk gas one is rig^d"^ ~ 0.0178. 

+00 of the two wall potentials are the same. In contrast, when ^w.lr = ^w,lj + o'w,lr, "wall 
potentials Vlr (dashed line) and Vlj (solid line) yield similar profiles except near the wall. 
These two last profiles can be made even closer to each other, including the area close to 
the wall, by changing the value of Zw.lr through ct^^lr (dotted line) and updating e„^LR 
accordingly. However, even if the profiles can become similar, important differences can still 
exist in the adsorption isotherms (cf next section). 

C. Isotherms 

Of particular interest is the dependence of the film thickness with respect to fi for a given 
T. This is crucial to understanding more involved configurations such as the contact angle 
case examined later. Because of the spatial variation of the density, it is in fact more relevant 
to consider an integral norm of the density, the "adsorption" , defined as 
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FIG. 3: Density profiles for /i < fisat at T = 0.7 Tc near an attractive wall. The wall potential 
is Vlj given by equation (jlOcp with parameters /JcCw = 1-8, = 1.25 d and = 0. Solid line, 
A* — /^sat = — O.OOI/cbTc; dotted line, /x — /Zgat = — 0.01A;bTc. In the former case, nn^d^ 0.684 which 
is close to the value of the density for z/d in the interval [3, 17]. 

where zi is the position of the dividing wall-fluid interface and rioo is in the present case equal 
to ngas (n = as 2 — — oo). The bifurcation diagrams for the isotherms are constructed 
with a continuation procedure which also involves an extra equation/geometric constraint 
which links the continuation parameter /i and n with the continuation step. Besides easily 
obtaining the density proflles for a wide range of values of /i, this approach has the additional 
beneflt that it allows the computation of unstable branches, which would have been more 
difficult if not impossible to obtain otherwise. 

Typical isotherms for the wall potential Vlj in equation (11 Pel) are depicted in figure [51 
For very attractive walls (for example, the case (S^e^ = 3 in figure [5]) and fi < /igat, the 
adsorption F is a monotonically increasing function of fi; in particular, only one solution 
exists for a given value of fi. At the saturation value /igat, the bulk liquid phase becomes as 
stable as the bulk gas one and the adsorption goes to infinity. For less attractive walls, a 
multi- valued S-shaped loop in the isotherm appears for values of fi say in [yU_,/i+] and with 
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FIG. 4: Density profiles for ^ — /Xgat — O.OOI/cbTc at T = 0.7 Tc near an attractive wall. Solid line: 
wall potential Vlj with /3cew = 1.8, = 1.25 d and = 0; dashed line: wall potential Flr with 
/JcCw = 1-8, cTw = 1.25 d and = o"„; dotted line: wall potential Vlr with /Sc^w ~ 3.52, a„ = d 
and = fw! dashed-dotted line: wall potential Vsr with /JcCw = 1-8, o"w = 1.25 d and = a^. 



three branches of solutions from which the middle one is always unstable (corresponding to 
variational maxima for This loop is associated with the presence of a first-order phase 
transition with respect to the adsorption: for a certain value of /i G three solutions 

exist, of which those in the lower and upper branches are equally stable (variational minima 
of VL of equal "depth"), and consequently the system can adopt for that value of /i two 
different adsorptions or equivalently film thicknesses. This transition is often referred to as 
the "prewetting transition" and a Maxwell construction in the variables (/i, F) can be carried 
out to find the coexistence value /icocx(r) of the chemical potential at which the two states 
have the same stability. For /i G /icocx], the lower branch is stable and the upper branch 
one metastable [VL has two variational minima but the one in the lower branch is "deeper" ) 
while for /i G [/icocx, /^+] the upper branch is stable and the lower one metastable. 

The value //sat imposes an upper bound on /icocx- We cannot have /icocx > /isat as a 
Maxwell construction in this region is not possible. As /icocx f^sat can have a Maxwell 
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construction with a thin film equally stable with an almost infinite thick film, while for 
= fJ'sat the most stable state is that of a thin film on the lower branch as is the case when 
the wall is only slightly attractive (e.g. (S^e^ = 0.8 in figure [S]). This is the signature of a 
partial wetting situation. Complete wetting on the other hand occurs whenever /icocx < /^sat 
and as /i — i- /igat on the upper branch. For a given value of e„, one can usually go from a 
partial wetting case to a complete wetting one by changing the temperature. 

Typical isotherms for the wall potentials Vsr and Vlr in equations f llOap and (llObp are 
shown in figures |6] and [TJ respectively. The isotherms appear qualitatively similar to those 
in figure |3 The effect of the potential ranges, however, becomes apparent in the upper half 
of the figures as, for a given isotherm, /i tends to figat as Td'^ goes to +00 much faster in the 
case of potential Vsr than in the case of potential Vlr. It may also be worthwile to note the 
important differences in the bottom half of the figures even though the distance from the 
wall is small there (i.e. the potential range is less critical): the values of /icocx and in 
the two cases differ significantly. This is confirmed in figure |8] where their dependance on the 
strength ew of the wall potentials is depicted. The intersection of /icocx with the x-axis (i.e. 
the /i = /isat axis) in figure [H] corresponds to the partial/complete wetting transition with 
respect to e^. For weakly attractive walls (e^ small), the wetting is partial and becomes 
complete for very attractive walls (ew large). In the present case, the transition occurs sooner 
for the short-range potential. 

Note that independently of the situation, complete wetting or partial wetting, there is 
always a film (albeit very thin) in contact with the wall. After all, the wall is always 
attractive. If it were repulsive, i.e. a hard or dry wall, there would be no contact between 
the fiuid and the wall. 



D. Square gradient approach 

A substantial simplification to equation f Ha|) can be achieved if we assume that n(r) varies 
smoothly in the range of 0p. Indeed, by expanding n(r) in a Taylor series and neglecting 
terms of 0(3) and higher, equation (Hal) becomes 

F[n] = F,[n] + j (^^\Vn{r)f - ^^^(r)^ (13a) 
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FIG. 5: Bifurcation diagram for F as a function of 6fj, = /U — fisat, for the attractive wall with 
potential Vlj (equation (jlOcp ). = and cjw = 1.25 d. The temperature is T = 0.7 Tc and the 
dividing interface is located at zj = z^. The wall parameter /ScCw varies from 0.8 (right curve) to 
3.0 (left curve) in steps of 0.1. 



where 

ir = -y y"drrVp(r). (13b) 

Instead of the integral equation ([7]) where F^\n\ is given by equation fjlbj) . we now have a 
differential equation: 

/ir(n) - Kl\n - Ian + ^(r) = /i. (14) 



This local approach has been widely used to compute ID density profiles (e.g. 2J]). It 
is also the starting point of Cahn-Hilliard- type/diffuse interface equations where the free 
energy is typically a function of density and its gradient (e.g. 13|). However, comparison 
with the DFT integral approach reveals a number of shortcomings, especially when it comes 
to interfaces. Figure IH] shows the liquid/gas interface obtained by solving equations ([7]) and 
(IT^ . Although the overall shape of the density profile is qualitatively similar for the two 
cases, differences appear in the transition area, where n(r) varies sharply. This is even more 
evident when a wall is present, as demonstrated in figure [TOl The local approach leads to 
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FIG. 6: Bifurcation diagram for F as a function of 6fj, = fi — fisat, for the attractive wall with 
potential Vsr (Eq. (jlOap ). = and = 1.25 d. The temperature is T = 0.7 Tc and the dividing 
interface is located at zj = z^. The wall parameter /ScCw varies from 1.2 (right curve) to 3.4 (left 
curve) in steps of 0.1. 

smoother curves with milder slopes and is unable to account for the small oscillations near 
the wall. 

Qualitative differences can appear between the two methods when considering the 
isotherms. For example, in figure [TT] the local approach shows a multi-valued curve and 
hence phase transition as opposed to a single-valued curve obtained from the integral ap- 
proach. 
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FIG. 7: Bifurcation diagram for F as a function of 6fj, = /U — fisat, for the attractive wall with 
potential Flr (equation (jlObp i. = and = 1.25 d. The temperature is T = 0.7 Tc and the 
dividing interface is located at zi = z^. The wall parameter f^^e^ varies from 1.2 (right curve) to 
3.4 (left curve) in steps of 0.1. 



IV. 2D PROBLEM 



A. Equations 



We now consider cases for which the translational invariance along one of the directions 
parallel to the substrate is broken: 



hoo 



+ ^ da;'dyn(a;', /)*2d \\j [x' - xf + {z' - zfj + V{x, z) = fi (15a) 
where 

/ + 00 
dy0p(/R^T7) (15b) 
-oo 

with 



R= ^J{x' - xf + {z' - zf. (15c) 

The x-dependence of the wall potential V could correspond, for example, to a chemically 
heterogeneous substrate. 
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FIG. 8: Special values (//_, ^coex and of the chemical potential as a function of e^. The solid 
lines correspond to wall potential Vsr and the dotted lines to wall potential Vlr. The remaining 
parameters are the same as in figures [6] and [71 

The integral in fllSb .c) can be performed analytically (see Appendix A). The compu- 
tations now are substantially more demanding compared to the previous ones, as besides 
being 2D, they also usually require a much larger number of iterations to achieve conver- 
gence. As in the ID case, we utilize a Newton scheme with an appropriately simplified 
Jacobian. Details are given in Appendix B. 

B. Prewetting transition 

As demonstrated in § IIII Ct for a range of values of temperature and chemical potential 
there exists a first-order phase transition with respect to the adsorption F, also known as 
the prewetting transition. The interface profile in the area joining two equilibrium films 
thicknesses at the prewetting transition is pictured in figure [12] for the long-range wall 
potential in equation f llOcp . The shape of the interface between the two films appears to be 
sufficiently smooth compared to the density profile of a liquid-gas interface as the interface 
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FIG. 9: Density profiles obtained by solving equations ([7]) (solid line) and (jl4p (dotted line) for 
T = 0.7 Tc and /i = ^isa.t{T). 

between the films is several tens of molecular diameters long, while in the latter case it is 
typically only of a few molecular diameters long. 

The shape of the transition area depends on the difference of the two co-existing equilib- 
rium thicknesses and also on the range of the wall potential. In figure [131 we present the 
profile obtained when the wall potential is a short-range one given in equation f llOal) . The 
transition now between the two films thicknesses appears to be more abrupt with a steep 
rim as the smaller thickness is approached giving the profile a pancake-type shape. 

C. Contact angle 

We now turn to the study of the density profile when /i = /isat and the wall is attractive 
with a potential given by equation f llOcp . Recall from our discussion in § IIII CI that up to 
A^coex ~^ fJ'sat "^^ cau have a Maxwell construction with a thin film co-existing with an almost 
infinite thick film, while at /i = /isat "we loose the Maxwell construction and the thin film is 
the most stable state. Again, as in § lIVBt two different conditions are imposed on each side 
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FIG. 10: Density profiles obtained by solving Eqs. ([7]) (solid line) and (jl4|) (dotted line) in the 
presence of a wall for T = 0.7Tc and fi — /.fsat(^) = — O-OOI/cb^c- The wall potential is Vi^j with 
parameters fSce^ = 1.8, = 1.25 d and = 0. 

of the domain. The wall parameter ew is chosen so that a stable film of finite thickness can 
be sustained on the wall. According to figure O this is true for Vlj when (3ce^ is at least 
in the range 0.8-1.4. The corresponding ID density profile computed in § IIIII will make up 
the boundary condition on the left side of the domain. On the right side, a thick film of 
arbitrary thickness /iright and shape is imposed (unlike the equilibrium film to the left, to 
the right we do not have an equilibrium film). For simplicity the latter is taken to be a step 
function: 



An equilibrium density profile for an attractive wall obtained under these conditions is 
depicted in figure [TH The main feature is the presence of a well formed angle between the 
two phases. This is more evident in the current profiles than in the profiles obtained in 
§ HVBl quite likely because now the film profile on the right side of the domain is no longer 
an equilibrium one. Note also the liquid film present in front of the contact line the 
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FIG. 11: Isotherms obtained by solving Eqs. ([7]) (solid line) and (|14p (dotted line) for T = 0.8 Tc. 
The wall potential is Vi,j with parameters /JcCw = 1-4, o"w = 1.25 d and z„ = 0. The dividing 
interface is located at zi = 



wall is attractive. Moreover, it is important to establish that the profile is not affected by 
the film thickness to the right. Indeed, as demonstrated in figure [13 the thickness of the 
film in the right side of the domain has little influence on the density profile since, for all 
values of /iright, the level curves are parallel to each other sufficiently far from the contact 
line. 

A cut in the density profile displayed in figure [14] far from the contact area and normal to 
the liquid- gas interface is shown in figure [TBI A comparison with the profile in figure [2] reveals 
that the two are very close to each other: the interface sufficiently far from the contact line 
is simply the usual liquid-gas interface, i.e. with no wall present. In other words, away 
from the contact line, the interface shape is not affected by the wall and it is like having the 
usual liquid-gas interface there. However, a deviation from the liquid-gas interface occurs 
in the three-phase area where a small incurvation of the density profile towards the wall is 
observed. We shall return to this point shortly. 
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FIG. 12: Equilibrium density profile between two equilibrium film thicknesses for the long-range 
wall potential Vlj given by equation (jlOcp . The temperature is T = 0.7 Tc and the chemical 
potential is = /^cocx- The wall parameters are /3cew ~ 1.5, = 1.2b d and = 0. The shade 
values provided to the right of the figure correspond to the level values of . 
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FIG. 13: Equilibrium density profile between two equilibrium film thicknesses for the short-range 
wall potential VsR given by equation (jlOap . The temperature is T = 0.7 Tc and the chemical 
potential is /i = ficocx- The wall parameters are /3cew = 2.0, fJw = 1.25 d and = 0. The shade 
values provided to the right of the figure correspond to the level values of nd^. 
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FIG. 14: Equilibrium density profile for a fluid in contact with an attractive wall, whose potential 
is Flj (equation (llOcp ). The temperature is T = 0.7Tc, and the chemical potential is ;U = /Usat- 
The wall parameters are P^e^ = 1.3, = 1.25 d and = 0. The contact angle is close to 37°. 



The surface tension of the interface is related to the grand potential Vt via 

^ — ^bulk 

where f2buik is the sum of the grand potentials of the two bulk phases and A is the area 
of the interface under consideration. Surface tensions for the planar wall/gas, wall/liquid 
and liquid/gas interfaces can be computed from the density profile constructed in figure UM 
Simple analytical expressions can also be derived in the sharp interface limit for an infinite 
system and are given in Appendix A. From the surface tension values, one can calculate the 
contact angle 6 from the Young equation, 

c^wg = o"wi + o-ig cos 6, (18) 

and contrast it with the one obtained by a direct geometric measurement in figure [17] using 
density level curves. We note that the choice of the specific level curve is not important since, 
as it has been demonstrate in figure [IHl away from the contact line the interface shape is 
very close to the one corresponding to the liquid-gas interface without wall, and consequently 
two different level curves would yield the same geometric contact angle. Results for several 
values of are given in Table [H A very good agreement is found between the two values of 
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FIG. 15: Density level curves (corresponding to n = (nnq + ngas)/2) for three different values of 
film thickness /iright imposed on the right side of the domain: height /d = 18 (right curve), 24 and 30 
(left curve). In all cases, the same film is imposed on the left side of the domain. All parameters 
are identical to those of figure [m 










^num ( ) 


^mcs ( 


1.0 


-0.0821 


-0.253 0.503 


70.2 


70.8 


1.1 


-0.106 


-0.360 0.503 


59.7 


60.7 


1.2 


-0.136 


-0.470 0.503 


48.4 


49.5 


1.3 


-0.171 


-0.582 0.503 


35.2 


36.6 



TABLE I: Wall/gas (cwg), wall/liquid (cJwi) and liquid/gas (dig) surface tensions, and computed 
(^num) and measured (^mcs) contact angles corresponding to figure [T71 



the contact angle. The main discrepancy between the two is quite likely caused by the error 
in the drawing of the dotted lines in figure [T71 The computations in figure [T7] also indicate 
that increasing the wall attraction decreases the contact angle (as expected). 
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FIG. 16: Liquid-gas interface (solid line) from Fig. [2] and cut in the density profile pictured in 
Fig. [HI normal to the liquid-gas interface (dotted line). The two profiles are very close to each 
other. 



D. Contact line 

A close inspection of the contact line area in figure [T7] reveals a deviation between the 
geometric profile suggested by Young's equation and the one obtained from DFT close to the 
contact line where the interface seems to bend towards the wall and away from the dotted 
lines corresponding to Young's equation. This feature appears to be in agreement with 



recent experiments by [22|. The experiments are actually done for a two-layer composite 
wall, but when the top layer is thin, the interface close to the contact line appears to be 
similar to that shown in figure [T71 



Level curves constructed with wall potentials Vg^ and Vlr, in equations f llOap and fllObp . 



are given in figures [18] and [191 respectively. The bend is not present in the first case which 
suggests that it is mainly due to the long-range nature of the attraction of the wall potential. 
A parallel can be made with the conclusions of § IIVBI regarding the density shape in the 
transition area between the two films: as the contact angle goes to zero, the profiles in 
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FIG. 17: Density level curves corresponding to n = {nnq + ngas)/2 for several values of e^. From 
right to left: /JcEw = 1-0, /Sc^w = l-l, /^c^w = 1-2 and /ScCw = 1-3. The contact angle in each case is 
measured by following the dotted lines. The wall potential is Vlj, the temperature T = 0.7 Tc and 
the chemical potential n = Usat- The wall parameters are = 1.25 d and = 0. 

Fig. [TS] seem to tend to the one in figure [12] giving rise to an abrupt transition between the 
two films while the ones in figure [12] seem to lead to the profile in figure [12] for which a 
smoother transition between the two films is observed. In contrast, no dramatic difference is 
observed between Vlj (figure [T7|1 and Vlr (figure [T9D implying that the shape of the repulsive 
part of the wall potential is less critical. 

Cuts in the contact line area normal to the wall are shown in figure [20] for the two 
potentials Vsr and Vlr in equations fllOap and fllObI) . respectively. We note that the density 
profiles are characterized by a small depression in the immediate vicinity of the wall {x/d ^ 
0.5), although this feature disappears in the case of the short range wall potential Vsr as we 
move away from the contact line towards the liquid. In contrast, it is present in the long 
range case even in the ID case as shown in figure [3] for 2; ?a 2d. 
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FIG. 18: Density level curves corresponding to n = (niiq + ngas)/2 for wall potential VsR and several 
values of £„. From right to left: /JcEw = 1-4, /ScCw = 1-5, /3ce„ = 1.6, fic^^ = 1-7 and /^cCw = 1-8. 
The temperature is T = 0.7 Tc and the chemical potential /i = /^sat while the wall parameters are 
(Tw = 1.25 d and = 0. 

V. SUMMARY 

We have examined the equihbrium of a fluid in contact with a sohd substrate through 
a DFT theory based on a perturbation approach in which the interaction potential is split 
into a repulsive and an attractive part. This leads to three distinct parts in the expression 
for the free-energy functional, corresponding to repulsion, attraction and external potential 
associated with the presence of the solid boundary, respectively. 

We first investigated the ID case of a liquid film in contact with the substrate. For values 
of the chemical potential less than but close to its saturation value, a liquid film exists in 
contact with the wall, even though in the absence of the wall the gas is the preferred state. 
The presence of the film is due to the attractive part of the wall potential. For a given 
temperature, we have constructed bifurcation diagrams for the film thickness as a function 
of the chemical potential. Such diagrams are typically S-shaped curves with three branches. 
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FIG. 19: Density level curves corresponding to n = (Tiiiq+ngas)/2 for wall potential Vi,^ and several 
values of e^- From right to left: /ScCw = 1-7, /Jc^w = 1-8, /Jc^w = l-9i f^c^w = 2.0 and /^cCw = 2.1. 
The temperature is T = 0.7 Tc and the chemical potential fi = /^sat while the wall parameters are 
(Tw = 1.25 d and = 0. 

Comparison of the DFT theory with a gradient approach in which the integral DFT equation 
is approximated with a differential one reveals some significant qualitative differences. Any 
gradient approach is by definition a local one and as such it neglects the non-local features 
induced by the intermolecular forces acting on the fluid system. 

The bifurcation diagrams reveal that there exists a special value of the chemical potential, 
which we referred to as "coexistence" value, in which a thin liquid film is equally stable with a 
thick one ( "prewetting transition" ). The interface that joins the two films is constructed with 
a fully 2D computation as now the translational invariance of the system in the direction 
parallel to the wall is broken. As the coexistence value tends to the saturation one, the 
thickness of the thick film tends to infinity. This then allows us to construct a liquid 
wedge in contact with the substrate and with a well-defined three-phase contact line. This 
wedge seems to persist for all distances from the substrate and hence it should eventually 
enter the macroscale. Comparison of the contact angle obtained from the DFT theory with 
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FIG. 20: Density profiles around the contact line area for the cases depicted in figure [18] (/3c = 1.8, 
solid lines) and in figure [19] (/JcCw = 2.0, dotted lines). Bottom/top curves correspond to the 
left /right of the contact line region in steps of 2 for x/d. 



a macroscopic model, such as the Young equation based on a mechanical force balance, 
shows good agreement a few molecular diameters far from the contact line. However, unlike 
the Young equation which requires information on macroscopic parameters, such as the 
surface tensions between the different phases, the DFT theory relies only on well-defined 
intermolecular parameters. 

Of particular interest would be the extension of this approach to: (a) Spatially hetero- 
geneous, chemical or topographical substrates. Such substrates have a significant effect on 
the wetting characteristics of the solid-liquid pair (e.g. 21|, [lot, jl9|, 20|); (b) The 
substantially more involved dynamic case. This would form the basis for the formulation of 
a dynamic contact angle theory from first principles without the need for any phenomeno- 
logical parameters. We shall examine these problems in a future study. 
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Appendix A: Expressions 
1. Surface tensions 

For a planar interface (i.e. the ID case), analytical expressions for the surface tensions 
can be derived in the sharp interface limit for an infinite system. They read: 

1 2 

cTig = --{nuq - ngas) / (Ala) 

(Twg = -2^gas^ + %as / dzV{z) (Alb) 







c^wi = -^nld^ + '^liq J dzV{z 



(Ale) 



where 



-oo 

1=1 dz I dz'^idiz' - z). (Aid) 

J-oo Jo 

Such expressions are also applicable in the contact-angle case sufficiently far from the contact 
line, i.e. when the interface is planar. 

In the following, we provide analytical expressions used both in the ID and 2D computations. 



2. Interaction potential 



The perturbative part of the interaction potential between fluid molecules in the Barker- 
Henderson case is: 

if r < cr 

. r / V r - 



Mr) 



4e 



if r > o". 



(A2) 
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a. Derivation of $ 



In the ID case, $, which enters the equation for the density, is defined as: 

$id(/3) = ^''°dxdy0p(v/x2 + i/2 + /32j . 

Let us introduce (3 = (3 /a. Using Eq. (lA2p . we obtain: 



Id 



1 /2 1 



if 1^1 < 1 
if 1^1 > 1. 



In the 2D case $ is defined as: 

$2d(/3) 

For < 1 we obtain: 



$2H = 8ea 



1 63 / TT 



256 1 2 



— arctan 



1-/32 



/32 



1-/32 /32(i28/3S + 144/36 ^ iQ8^4 ^ 210/32 + 315) 



1280 



/35 ] 8 2 



arctan 



'1-/32 



/32 



1-/32 /32(2/32 + 3) 
~^2 § 



and for |/3| > 1: 



3 1 

2 /35 



21 1 
32^ 



(A3) 



(A4) 



(A5) 



(A6) 



(AT) 



3. External potential 

The Lennard- Jones expression of the wall potential can be obtained by considering that 
the wall is made of particles which interact with the fiuid molecules via a pair potential of 
the Lennard- Jones form: 

^■m-<((t)"-(t)°) 

where r is the distance between a wall particle and a fiuid particle. 

Let us consider the case of a planar wall. Our coordinate system has the three unit 
vectors (i, j, k) with k an outward unit vector normal to the wall. The wall is described as a 
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continuous medium and the wall particles density, denoted by n^, is taken uniform. The the 
edge of the wall, which is in direct contact with the fluid, is located at z'^ so that ra^ = if 
z > z'^. The overall potential exerted by the wall at position (xq, 2/0, -^o) in the fluid domain 

(^0 > O is, 



V[xQ,yQ,ZQ) 



dr n^v{r) 



(A9) 



^<Zw 



where rzw denotes the wall particles density and is assumed constant. As the wall is invariant 
in the (i, j) directions, V only depends on zq. Using cylindrical coordinates, we obtain. 



2n 



V{zo) 

which for zq > z'^ yields: 



RdR / d^ / dzn„v[ B? + (^0 - z)' 
Jo 



Zq — 



15 \ 2:0 — z'^ 
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Appendix B: Numerical method 



1. ID case 



The equation to solve reads: 

/+00 
dz'n{z')^iA{z' - z) + V{z)= (Bl) 
-00 

for n{z) in the domain =\zq^zn^[ with conditions n = ioi z < zq and n = for 
z > zj\f^. In the presence of the wall we assume Zq > z„ (or zq > z'^). We introduce a uniform 
mesh: Zi = Zo + Az i, i = O...Nz. The integral is computed by using a trapezoidal rule in the 
interval [zq, zn^] while analytical expressions are utilized outside this interval, obtained by 
using (lA4p and the fact that the density there is constant (= n\^g)- We then obtain a system 
of Nz — 1 nonlinear equations for the unknowns rii = n{zi), i = l...Nz — 1, which are solved 
by using a Newton scheme. In order to speed up the scheme, the off-diagonal elements of 
the Jacobian matrix are neglected as they are smaller compared to the diagonal ones (this is 
especially so far from the diagonal). The Jacobian matrix can then be inverted analytically 
in each iteration step. 
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2. 2D case 



The equation to solve reads: 

/ + 00 
dx'dz'n{x', z')<^2d (V(^'-xOV(^'-^)^) + V{x, z) = fi (B2) 
-oo 

for {x,z) G 'D2d where =]xo,X]\f^[x]zQ, zj\f^[. In the presence of a wall, Zq > z^ (or 
zo > ^w)- The conditions outside this domain are n{x, z) = ni{z) if x < xq, n{x, z) = Uj^lz) 
if X > xtv^, n(x, z) = Uf, a z < zq and n{x, z) = Ut ii z > zjsi^ The functions ni and rij. are 
such that ni{z) = Uf, and = rib ii z < Zq, and = ra^ and nj.{z) = rit ii z > Zi\f^. 

This choice of conditions for x < xq and x > xat^ allows us to examine, for example, the 
case of a fluid in contact with two liquid films of different thicknesses, one on each side. 

A uniform mesh is used to compute n(x, z): Xj = Xq + Ax j, j = O.-.N^, and Zi = z^ + Azi, 
i = O...Nz. Equation ( 1B2I) is solved for riji = n{xj,Zi), j = l.-.N^ — 1 and i = l...Nz — 1. 
The integral appearing in the equation is computed by dividing the integration domain into 
three parts: 'D2d, 1^2d ~ ^2d and ] — oo, 4-oo[x] — oo, +oo[— Pgd- The domain "Dgd is defined 
as =]3^0'^7V^[^]4' -^TvJ' where Xq = Xq-ct, x^^ = Xn^ +a, z'q = ZQ-a and z'j^^ = zn^+ct. 
Outside Pgd' Eq.(jX7]) is used for any riji G 'D2d and analytical expressions can be derived 
for the integral when |x', z'\ — > oo where the density is constant (= ni^g). In Pgd ^ ^sd, the 
integral is computed numerically by using either Eq. ( ]A6[) or Eq. ( ]A7P once for all at the 
start of a run. In V2d, the integral has to be carried out numerically at each iteration. 

The discretization of Eq. (\B2\i leads to a set of {N^^ — 1) x (A^^ — 1) nonlinear equations 
for the [Nx — 1) x (A^^ — 1) unknowns nji. It is solved by using a Newton scheme. To 
speed up the computations, the Jacobian is made sparse by neglecting contributions when 
|ri — r2|, where ri 2 the positions of two particles, is larger than a few molecular diameters: 
the non-diagonal terms of the Jacobian which involve $2d(ri — r2) (for the diagonal terms 
ri — r2 = 0) are getting small as we move further from the diagonal since |ri — r2| increases 
so that $2d decreases. By doing so, more iterations are needed for convergence but each 
one is considerably faster than inverting the full Jacobian. Finally, at each iteration step, 
we have to solve a linear system involving a sparse matrix. This stage is performed with a 
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Gauss-Siedel method. 
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